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Abstract — In a distributed space-time coding scheme, based on 
the relay channel model, the relay nodes co-operate to linearly 
process the transmitted signal from the source and forward them 
to the destination such that the signal at the destination appears 
as a space time block code. Recently, a code design criteria for 
achieving full diversity in a partially-coherent environment have 
been proposed along with codes based on differential encoding 
and decoding techniques. For such a set up, in this paper, a non- 
differential encoding technique and construction of distributed 
space time block codes from unitary matrix groups at the source 
and a set of diagonal unitary matrices for the relays are proposed. 
It is shown that, the performance of our scheme is independent 
of the choice of unitary matrices at the relays. When the group 

. is cyclic, a necessary and sufficient condition on the generator 
of the cyclic group to achieve full diversity and to minimize 
the pairwise error probability is proved. Various choices on the 
generator of cyclic group to reduce the ML decoding complexity 
at the destination is presented. It is also shown that, at the source, 
if non-cyclic abelian unitary matrix groups are used, then full- 
diversity can not be obtained. The presented scheme is also robust 

\ to failure of any subset of relay nodes. 

Index Terms — Cooperative communication, cyclic groups, dis- 
tributed space-time codes and unitary space-time codes. 

I. Introduction and preliminaries 

Co-operative communication can be based on a relay chan- 
nel model where a set of distributed antennas belonging to 
■ multiple users is exploited for achieving spatial diversity JT|- 
|H). In 0, the idea of Space-Time Coding (STC) devised for 
point to point co-located multiple antenna systems is applied to 
a two-hop wireless relay network and is referred as Distributed 
Space-Time Coding (DSTC). The technique involves a two 
phase protocol where, in the first phase, the source broadcasts 
the information to the relays and in the second phase, relays 
linearly process the received signals and forward them to the 
destination such that the signal at the destination appears as a 
Space-Time Block Code (STBC). 

In the above technique, the destination may or may not have 
the knowledge of channel fade coefficients (i) from the source 
to the relays and (ii) from the relays to the destination when de- 
coding for the source signal. DSTC when the destination does 
not have the knowledge of the channels from the source to the 
relays but has the knowledge of the channels from the relays 
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to itself is called partially-coherent DSTC, which is the subject 
matter of this paper. For a partially-coherent set up, in Q, a 
code design criteria is proposed for achieving full diversity. A 
differential encoding and decoding strategy with a class of full 
diversity achieving codes using cyclic division algebras is also 
proposed. Inspired by non-coherent MIMO unitary differential 
modulation [6|, differential distributed space time codes have 
been proposed in |7|, [8] and J9)- 

In J5), it is shown that the coding problem for a partially 
coherent setup is to distributively design a finite set of unitary 
matrices such that the matrix obtained from juxtaposing any 
two distinct matrices from the set must be full rank. In other 
words, the problem is to design non-intersecting subspaces 
such that the principal angles between subspaces is as large 
as possible. However, with the use of differential techniques 
flU, 10, El, ID, the coding problem gets transformed in to a 
problem of constructing a finite set of unitary matrices such 
that the difference matrix of any two matrices is full rank. This 
is a well known design criteria for relay networks with amplify 
and forward protocol in a coherent detection environment and 
code constructions based on the above criteria are available in 
the literature. In this paper, we propose a method to explicitly 
construct non-intersecting subspaces for a partially-coherent 
set up. i.e we construct a set of unitary matrices in a distributed 
way such that the matrix obtained from juxtaposing any two 
distinct matrices from the set must be full rank. 

The contributions of this paper can be summarized as 
follows: (i) for the partially-coherent set up, we introduce 
a non-differential coding scheme and obtain an expression 
for ML decoding metric (Theorem [TJ for the general class 
of unitary Distributed Space-Time Block Codes (DSTBCs) 
(See Definition |2). (ii) We construct unitary DSTBCs using 
cyclic unitary matrix groups called Non-differential Cyclic 
Distributed Space-Time Codes (NCDSTC) (Definition 0J and 
provide a necessary and sufficient condition on its generator 
such that the NCDSTC is fully diverse (Theorem 0I and 
the pairwise error probability (PEP) is minimized, (iii) Also, 
we provide conditions on the choice of the generator of the 
cyclic group so as to reduce the decoding complexity at the 
destination, (iv) We show that the proposed scheme is robust 
to failure of a subset of the relay nodes. 

The remaining part of the paper is organized as follows: 
In Section [II] along with the signal model, partially-coherent 
distributed space time coding technique is briefly reviewed 
and an ML decoding metric for a class of unitary DSTBCs 
is presented. The scheme of non-differential unitary DSTC 
from cyclic unitary matrix groups is introduced in Section 
ITTT1 where we provide a necessary and sufficient condition on 
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Fig. 1. Wireless relay network model 



the generator of the cyclic group to achieve full diversity 
and to minimize the PEP. In Section |IV| an ML decoding 
metric for the proposed codes is presented along with two 
reduced decoding complexity ML decoders. Details on the 
robustness of our scheme to the failure of a subset of the relay 
nodes is also presented. Possible directions of future work and 
concluding remarks constitute Section [V] 
Notations: For a complex matrix X, the matrices X*, X T , 
X H , |X|, Re X and Im X denote, respectively, the conjugate, 
transpose, conjugate transpose, determinant, real part and 
imaginary part of X. For a complex matrix Y of the order same 
as X, X © Y denotes the Hadamard product of X and Y. The 
T xT identity matrix is denoted by It and Ot denotes the T- 
length vector of zeros. We use \x\ to denote the absolute value 
of the complex number x and E [x] to denote the expectation 
of the random variable x. We write x ~ CSCQ (fi, T) when x is 
a circularly symmetric complex Gaussian random vector with 
mean \i and covariance matrix T and use j for The set 

of integers and the set of complex numbers are, respectively, 
denoted by Z and C. 

II. Partially Coherent Distributed space time 

CODING 

A. Signal model 

The wireless network considered in Figure Q] consists of 
R + 2 nodes each having single antenna which are placed 
randomly and independently according to some distribution. 
There is one source node and one destination node. All the 
other R nodes work as relays. We denote the channel from the 
source node to the j-th relay as fj and the channel from the j- 
th relay to the destination node as gj for j = 1, 2, • ■ • , R, The 
following assumptions are made in our system model: (i) all 
the nodes are subjected to half duplex constraint, (ii) fading 
coefficients fj,gj are i.i.d CSCQ (0,1) with coherence time 
interval, T. (iii) all the nodes are synchronized at the symbol 
level, (iv) Destination knows the fading coefficients gj's but 
not fj's. 

Every transmission from the source to the destination com- 
prises of two phases. In the first phase the source trans- 
mits a T length complex vector from the codebook S = 



{si, S2, S3, • • • ,Sl} consisting of information vectors S; G C T 
such that E [sfs;] = 1, so that P\T is the average transmit 
power. In particular, P\ is the average transmit power used at 
the source node for every channel use. When the information 
vector s is transmitted, the received vector at the j-th relay 
is given by Yj = y / P^Tf j s + rij, j = 1, 2, ■ ■ ■ , R. In the 
second phase, all the relay nodes are scheduled to transmit T 
length vectors to the destination simultaneously. Each relay 
is equipped with a fixed T x T unitary matrix Aj and is 
allowed to linearly process the received vector. The j th relay 

rA,r, where P2 is the 



is scheduled to transmit tj = y ^ 1+Pi ^ j.j 
average transmit power used at each relay for every channel 
use. The vector received at the destination is given by 



R 



' P1P2T 
(1 + P1) 



Sh + n 



(1) 



where w ~ CSCQ (0,It) is the additive noise at the desti- 



nation, n = 



P-2 



JT+FY) ' w The equivalent 

channel h is given by h = [f 1 g 1 } 2 g 2 • ■ ■ /rSrF € C R . The 
codeword matrix S is given by S = [AjS A 2 s ■ • • A^s] € 

C TxR . 

Definition 1: The collection C = {[AiS A2S ■ • ■ A_rs]} 
of codeword matrices when s runs over S, is called the 
Distributed Space-Time Block code (DSTBC). 



B. ML Decoder 

Since the relay matrices Aj are unitary, the random vectors 
w and iij are independent and Gaussian and since gj are 
known at the receiver, n is a Gaussian random vector with 

E [n] = Ot and 

Assume that S € C. When the destination has the knowledge 
of gj's but not of fj's, y is a Gaussian random vector with 

E[y\S,gj]=0 T and 
E [yy^lS, gj ] = S y = P SS,S ff + 7 I r 



where, 



' 7=(l 1 + I TTk)Ef = i(l5,| 2 )) and 
. S ft = diag(| 5l | 2 ,---| fffl | 2 ). 
Then, the conditional pdf P(y|S, gj) is given by 

P(y\s,g,) = ^ex P (-(y H i:- 1 y)). 

The partially-coherent ML decoder decodes to a codeword S 
where 

S = arg max. P(y\S, gj). (2) 

Definition 2: A DSTBC is called a unitary DSTBC if 
every S G C in Definition Q] satisfies the condition that S H S = 
Hr, where t is a constant real number independent of the 
codeword S. 
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Theorem 1: For a unitary DSTBC, the partially-coherent 
ML decoding is given by 

arg max (y ff SGS H y) , (3) 

where G = diag • • • (i R ) and 

^ = l%| 2 (lft| 2 +7^ 1 )" 1 for j = l,- ,R- 

Proof: Using the well known result, |I + AB| = |I + BA| 
and the appropriate use of the matrix inversion lemma in (0, 
the result follows. ■ 



C. Design criteria 

Chernoff bound on the PEP for the decoder in (f2| is given in 
|5| following which a criteria for designing unitary DSTBCs 
in order to minimize the PEP is also provided. For achieving a 
diversity order of R, the term |S nm S„ TO | has to be non zero for 
all to, n such that n ^ to, where S nm € C Tx2R is obtained 
by juxtaposing two codewords S„ and S TO as, 



(4) 



Further, for large values of P (where P = Pi + RP2 is the 
total power used at all nodes for every channel use), PEP can 
be minimised by maximising |S^ m S Ilm | for all n 7^ to. A 
necessary condition for this is T > 2R. Hence, in the rest of 
the paper, we consider T = 2R. 

In the following section, we construct codes that satisfy 
the above design criterion for full-diversity using Generalized 
Butson-Hadamard matrices defined as follows: 

Definition 3: A Generalized Butson-Hadamard (GBH) ma- 
trix ifTOl is a T x T matrix M with entries such that 

MM fl = M H M = TJ T 

and the conjugate of every entry rriij of M is its inverse to" 1 

* -1 
i.e TOjj = rrij • . 



III. A NON-DIFFERENTIAL SCHEME AND CODES USING 
CYCLIC UNITARY MATRIX GROUPS 

We introduce a non-differential encoding scheme using a 
cyclic unitary matrix group [;6| of diagonal matrices at the 
source and relay matrices constructed from GBH matrices j8|, 
IfTOl . Appropriate ingredients required at the source and the 
relays are as follows: 



A. At the Source 

The source node is equipped with a finite abelian group U 
of size L consisting of 2R x 2R diagonal unitary matrices and 
a 2R x 1 complex vector, x. Let the group U be written as a 
direct product of K cyclic groups U v , each of order L v for 

v = 1, 2, • • • , K as, 



U = U\ x U2 x • • • Uk 



(5) 



where L — JlfLi L v . When the set of L„'s for v = 
1,2, •■■ ,K are pairwise relatively prime, the group U will 



be cyclic. Using the mixed-radix notation I = (li, I2, ■ ■ ■ , Ik), 
where each 1 < l v < L u , we denote every element of U as 

= n ^ 

where D„ is a generator of the v th cyclic group U v given by 

(6) 

For each v, u lv e {1, 2, • • • L v — 1} for i = 1, 2, • • • , 2R 
have to be chosen appropriately to achieve full-diversity and 
to maximise coding gain. The i-th diagonal element of D' is 
given by 



cxp 



Hence, the source is equipped with a code {U,x} where x = 
1 1 ••• 1] T . 

V2R 1 1 



B. At the relays 

Let M be the Rx R GBH matrix, usin^ which we construct 
a 2R x R matrix T as T = [M T M t ] . The columns of T 
are used to construct diagonal unitary matrices for the relays 
as Aj = diag (Ty, r 2j , • • • T 2 rj) where T tJ is the (ij)-th 
element of T and j = 1, 2, • • • , R. 

C. Encoding Scheme 

The source node maps log 2 L bits of information on to one 
of the L matrices from U say, D fc for some k where 1 ^ k ^ L 
and transmits the 2R x 1 vector D fe x to all the relays. Each 
relay performs linear processing on its received vector using 
the unitary matrix Aj and transmits a 2R length vectors to 
the destination. The distributed space-time codeword is of the 
form, 



[A X D* 



A 2 D fc x 



D [Aix A 2 x 
1 -»fc n 



■ • Ai?D fe x] 
A^x] 



(8) 



'2R 



:D r 



where we have used the fact that D and A^ commute for all 
j. It can be observed that S^S^ = Ir. 

Definition 4: The collection C of 2RxR codeword matrices 
shown below when k runs from 1 to L 

1 



C 



'2R 



=D fe r 



(9) 



is called a Non-differential Abelian DSTBC (NADSTBC) and 
if the group is cyclic, then it is called a Non-differential Cyclic 
Distributed Space-Time Block Code (NCDSTBC). 

D. Design criteria for D u 

In this subsection, we derive a criteria for choosing the 
generators of each of U v such that a NADSTBC is fully 
diverse. 
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( ^ l v {u gv - \ ( 1 V {U ZV -U( R+g ) v ) 

exp \j2tt I + exp I j2tt ^ j-± — 



(7) 



Theorem 2: The diversity order of a wireless relay network 
with R relays using a NADSTBC is R if and only if the gen- 
erators Y) u in © are chosen such that, for each i = 1, 2, • • • R, 



K 

E 



Z 



(10) 



for all pair (1, 1) such that I ^ I. 

Proof: Let -J=D Z I\ 
as in Q^p is given by, 

S,r = - 



7=D'r e C. The 2R x 2i? matrix 



/2i? 



DT DT 



Since ISffS 



irwjji — l^n'S^I, we consider the choice of D„ for 
v=l,2,---K such that |S ;; -Sf | ± where 



1 

2~R 



Vn Via 
V 21 V 22 



(11) 



with Vn, V12, V 2 i and V22 are Rx R diagonal matrices given 
by 



Vn — 2R I RxR ; V 2 i=Vf2; V 2 2 — 2R I Rx r 

and for each z from 1, 2, • ■ • , R the z-th diagonal entry of V12 
is given in ©. Applying the result, 



A B 
C D 



IAIID - CA : B 



where A,B,C and D are square matrices of same order and 
A -1 exists, on the matrix in (fTTT i and using (O, we get 

2R 



ISftSfl 



1 
2R 



|V 11 ||V 2 2-V2iVr 1 1 V 12 



1 



= ' 2R ' ' diag ^ 62 



Sr)\ 



where, 



R 



2 ^E 



For |S H -S„- 1 not to be zero, the product <5i#2 ■ • • 6 R must be a 
nonzero value. For each Si not to be zero, tij„ and u^ +R ) v 
for all v = 1, 2, • ■ ■ , K has to be chosen such that 



K 

E 



(Zy — lu){Ui u — U( R+i ) v ) 



for all 1,1 such that I 7^ I. It is clear that the above condition 
on D„ is both sufficient and necessary. ■ 
Corollary 1: DSTBCs from non-cyclic abelian groups are 
not fully diverse for any choice of D„. 

Proof: Since each of the cyclic groups U u C U, from the 
special case of Theorem [2] with K = 1, every Dj, must satisfy 



the criteria, 



gcd (Uui - u u u 



(i+R),L v ) 



1 for i = l. 



R. 



■ ■ ■ exp 



Since at least two of the L v share a factor and each D„ satisfy 
the above criteria, condition in dTOb of Theorem [2] will not 
be satisfied for at least one pair of values I and I. Hence, 
non-cyclic abelian DSTBCs cannot provide full diversity when 
used in the proposed scheme. ■ 

Henceforth, throughout the paper, we will consider only 
NCDSTBCs. We explicitly state Theorem for NCDSTBC 
as a corollary as: 

Corollary 2: Let U in (|5]l be a cyclic group with a generator 
D given by 

L 

A necessary and sufficient condition on the choice of D such 
the that diversity order of the wireless relay network is R, is 
given by 

gcd(ui — u R+ i, L) = 1 for each i = 1, 2, • • • R. 

Proof: The result follows as a special case of Theorem |2] 
for K =1. ■ 

Specializing the proof of Theorem |2] for the cyclic group 
case, the coding gain for NCDSTBC is given by 

R 

1 — cos 



D = diag <^exp ^j'27r— J , exp \J2tt- 



( -o U2R \ 
{ j27r —) 



} 



min 



2tt(1 - l)(ui - u i+R ) 



L 



(12) 

Thus, the vector v has to be chosen such that <fr is maximized, 
where 



(tlx - Ur+1,U2 - u R+2 •••Ur— u 2 r) . 



(13) 



From the result of Theorem[2]and Corollary [2] the performance 
of NCDSTBC is characterized by ( TT2l which is independent 
of the unitary matrices Aj at the relays and dependent only 
on the generator of the cyclic group. 



E. On the choice of T 

Instead of a GBH matrix, M can also be any Rx R unitary 
matrix with no zero entry, using which T is constructed as 

" M " 
M 

Though, the relay matrices Aj are not unitary in this case, 
simulation results show that, for a given U at the source, BLER 
(Block error rate, which corresponds to errors in decoding 
a codeword is considered as error events of our interest 
throughout the paper) performance is independent of the above 
choice of M. This can be seen from Figure |2] which is obtained 
using the code parameters L = 16 and v = [5, 7, 11, 1] for a 
network with 4 relays. However, decoding metric in (O doesn't 
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- With GBH matrices 

- With unitary matrices with no zero entry . 




Fig. 2. Plot illustrating the independence of the BLER on relay specific 
diagonal matrices 



Fig. 3. BLER plot for 4 relays with L = 16 for different v 



hold in the latter case. Therefore, we continue to use only GBH 
matrices to construct Aj in the rest of the paper. 

IV. Reduced ML Decoding complexity for 
NCDSTBCs and Node failures 

For a general class of unitary DSTBCs used in a partially- 
coherent scheme, an ML decoding metric is given in (O. In 
this section, we specialize this decoding metric for the class 
of NCDSTBCs. 

Theorem 3: For a NCDSTBC in ©, the partially-coherent 
ML decoding given in (0) reduces to 



arg fc max ^ (yf + Re {yf G 2 y 2 } + yf G 3 y 2 ) 



(14) 



where, 

y± = [yi V2 

d = Ti O, G 2 = T 2 ft and G 3 = T 3 ft 



yn] T , y 2 = [yn+1 vr+2 ■ ■ ■ y2R] T and 



where denotes the Hadamard product and Ti,T 2 , T3 € 

£ RxR giyen by [ Ti ]_ = w ( Ui -u s )k^ [r 2 ] i 3 = 

[T 3 ] i • = w ( - UR +-- UR +i )k and ft is the R x R Hermitian 
matrix, with the th element given by 



ft 



1,3 



(15) 



l} are the entries of M, the chosen GBH matrix and 
9i\ 2 (\9j\ 3 +'TP-T 1 ^ j = !,-■■ ,R. 



where m 

Pi^WW+lp- 1 ) ' for 3 = 1 

Proof: The result follows by substituting (O in d3) 
In the rest of this section, we present various choices on 
the vector v in <fl~3T > such that the computations required in 
decoding using the metric given by ( fl4l i is reduced. For large 
value of P, performance of the scheme is determined by the 
choice of the vector v. There are several ways of choosing 
the vector v from various choices of u = [iti, it 2 • • • u 2 r]- The 
optimal choice of v can still be made with appropriate values 
of u\ , U2 ■ ■ ■ ur by keeping 



From the above choice of v, T 3 will be independent of the 
exponent k. Therefore the ML decoding metric in (IT4~b reduces 
to 



arq max 

k 6 1,2,— L 



(ylGiy 1 +Re{y t 1 G 2 y 2 }). (16) 



Hence, compared to ( [Pil l, the number of computations required 
in ( fl6b to decode a codeword is reduced. 

The vector v = (ui,f 2 , • • • ,vr) can also be chosen to be 
of the form 

v\ = v 2 = ■ ■ ■ = v R , 

by choosing u\ = U2 = ■ ■ • = ur and u_r+i = w/?+ 2 = 
■ • ■ = U2R such that gcd (u\ — Uij+i, L) = 1. In such a case, 
Ti and T 3 are independent of the exponent k and hence the 
decoding metric is written as 



arq max 

k g 1,2,- 



Re 



rfny 2 } 



(17) 



In general, in order to estimate the most likely transmitted 
codeword from the source, the destination needs to perform 
matrix multiplication of order R x R according to ( fl4l ) for 
every k from 1 to L, where as, in the case of the metric 
given in ( TTTb . the destination needs to perform complex 
number multiplication for every k to estimate the most likely 
transmitted codeword. Hence, the number of computations for 
decoding a codeword at the destination reduces substantially. 

Figure [3] shows the BLER performance against the total 
power, P for different choices of generator D with L = 
16 for a network with 4 relays. The plot shows that, for a 
choice of v with vi = V2 = ■ • • = vr, the NCDSTC loses 
out close to 5-6 db compared to the one with the choice of 
v = [11, 7, 3, 1]. The relay matrices used in all the simulations 
are Ai = I 8 , A 2 = diag {1, -1, 1, -1, 1, -1, 1, -1}, 
A 3 = diag{l, 1,-1, -1,1, 1,-1,-1}, andA 4 = 
diag {1, — 1, — 1, 1, 1, — 1, — 1, 1}. The decoding metric 
given in (O has been used to carry out the simulations. In all 
our simulations, the total power P has been distributed to the 



source node and the relay nodes as Pi 



and Po 



p 

2R/ 



UR+1 — Ur +2 — Ur +3 



U2R- 



A. Failure of relay nodes 

In this subsection, we analyze the behavior of the proposed 
codes when a subset of the relay nodes fail to co-operate with 
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techniques such that the unitary distributed space time block 
codes can be constructed algebraically. 



Fig. 4. BLER performance of NCDSTC with one relay down 



the source i.e when a subset of the relays do not transmit 
information of the source to the destination. The structure of 
every codeword Sk from the NCDSTC C is given in Definition 
|U Similar to codes designed for non-differential non-coherent 
collocated MIMO systems, the columns of a codeword Sk 
can be viewed as a basis of an i?-dimensional subspace 
of C 2R . Two codewords satisfy the full diversity condition 
means that, the i?-dimensional subspace spanned by columns 
of each of them are non-intersecting. Therefore, a fully diverse 
NCDSTC can be viewed as a finite set of non-intersecting R- 
dimensional subspaces in C 2R . From the structure of the code, 
it is easy to see that, every relay contributes a basis element to 
each of the subspace (code-word). Let denote a codeword, 
after say, d of the R relays fail to co-operate with the source. 
The columns of a new codeword Sfc span a (R—d) dimensional 
subspace of _R-dimensional space spanned by the columns of 
Sfc. Hence, the new set of subspaces (codewords) remains to 
be non-intersecting and the new code C will remain to be fully 
diverse with diversity gain R — d. Figure [4] shows simulation 
results on the behavior of BLER when one of the relay fails 
in a network with 4 relays. The plot is obtained by using a 
code with parameters L = 16 and v = [11, 11, 11, 11]. 
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V. Discussion 

We considered the problem of code construction based 
on non-differential technique for wireless relay networks in 
a partially-coherent environment. Unitary distributed space 
time block codes were constructed using a cyclic group of 
diagonal matrices at the source and diagonal relay matrices 
constructed from GBH matrices. We have shown that the 
BLER performance of the proposed scheme is independent 
of the relay specific unitary matrices. Further, a necessary and 
sufficient condition on the choice of the generator of cyclic 
group is provided for achieving full diversity and to minimize 
the PER Various choices on the generator of cyclic group to 
reduce the ML decoding complexity at the destination was 
provided. The resistance of the proposed scheme to the failure 
of a subset of the relay nodes is also verified. It was also shown 
that, using non-cyclic abelian group of diagonal matrices 
at the source doesn't provide full diversity. An interesting 
direction for future work is to design non-differential encoding 



